Quasiparticle Excitation Spectrum of an Isolated 
Vortex in a High— Temperature Superconductor 

Irena Knezevic and Zoran Radovic 
Department of Physics, University of Belgrade, 
P. 0. Box 368, 11001 Belgrade, Yugoslavia 



Abstract 

The quasiparticle energy spectrum of an isolated vortex in a clean layered <i-wave su- 
perconductor is calculated. The Bogoliubov-de Gennes equations are solved perturbatively, 
within the model of step- variation of the gap function, adjusting the Caroli, de Gennes and 
Matricon approach for low-lying excitations in cuprates. A large peak in the density of states 
in a "pancake" vortex is found, as a consequence of two-dimensionality and strong coupling. 
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1 Introduction 

Enhancing interest has been focused on quasiparticle spectra properties around vortices. 
Density of localized states near the Fermi surface in the vortex core in conventional super- 
conductors (s-wave pairing, weak coupling, isotropic) is equivalent to the electron density of 
states of a normal- metal area with radius of the order of the coherence length £ ~ 10 nm.[l- 
3] High-temperature superconductors (HTS) are layered, having a cylindrical Fermi surface, 
the superconductivity is of <i-wave and strong coupling type, the coherence length being 
much smaller, £ ~ lnm. This implies a crucial difference in the quasiparticle spectra in the 
vortex core between conventional superconductors and HTS. 

Numerous computations of the quasiparticle spectra around two-dimensional (2D) vor- 
tices have been performed recently. [4-10] The common feature of layered superconductors is 
the appearance of a large peak in the quasiparticle density of states (DOS) in the vortex 



core, in agreement with the scanning tunneling microscopy (STM) observations on NbSe 2 []TT 



and YBCO.[12] However, in BSCCO the gap structure is found, without any signature of 



the quasiparticle bound states in the vortex core. [|13 | 

In this paper, the Bogoliubov-de Gennes equations for the vortex core are solved 
analytically, using the perturbation theory within the model of step-variation of the gap 
function. Since the bound states spectrum is not quasicontinuous for 2D and strong coupling 
superconductivity, this approach modifies the Caroli, de Gennes and Matricon results for 
low-lying states. Our results may confirm that DOS has either one large maximum or gap 
structure, depending on the value of k F ^ in HTS. 
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2 The Model Solution 



For an isolated "pancake" vortex in HTS in a low magnetic field applied along the c-axis, 
using the cylindrical coordinates (r,ip,z), z \\ c, we take the gap function 



A = A(r,0)e~^, 



(1) 



due to the flux quantization, neglecting the vector potential contribution, A ss 0. 

Furthermore, neglecting the self-consistency condition, we assume the spatial variation 
of the gap function 



A = A(r,9) = 



0, r < r c (6) 
A(0), r>r c (9) 



(2) 



where 



A(0) = A cos(2#) 



(3) 



for d-wave pairing, 9 being the polar angle in k-space. The vortex core radius r c {9) is of the 
order of the BCS coherence length £o(#) — hv F /n\A(9)\. 
Taking for the Bogoliubov amplitudes u and v 

(u(t,0)\ ( f + (r,9)e< v - 1 ^' 



\v(r,9) 



\f-(r,0)e 



i{v+l/2)<p 



(4) 



where 2u is an odd integer, we write the Bogoliubov-de Gennes equations in the form 



(H + H 1 )f = Ef, 



n = -La z + A(r, 9)a x , H x = 



uh 2 
2mr 2 ' 



/ 



u 



(5) 
(6) 



where ex,- are the Pauli matrices. 



We solve Eqs.(Bf) perturbatively. First, [EUJ the eigenproblem Hof = E f, i.e 



0. (8) 



(L + E -A(r,6)\ / / 4 
K A(r,9) L-E ) V/_ 



Second, the energy shift according to the first-order perturbation theory. 

We restrict our consideration to localized states, with energies below the bulk gap, 
I-EqI < |A(0)|. The solutions of Eq. ([8]) are real-valued, finite for r = 0, and zero for r — > oo, 

f± = A ± J n (q ± r) , (9) 

for r < r c (8), and for r > r c (8) 

f + = BHV(kr) + B*HW(k*r), (10) 



/- = -* yA2 ^ g) E ° (BH*\kr) - B*H%\k*r)) + (BH^(kr) + B*H^\k*i 

(11) 

Here, J and if are the Bessel and the Henkel function, respectively, A± are real and B is a 



complex constant, q± = \J k\ ± 2mE /h 2 , and = \jk 2 F + ty2m(A 2 (6) — Eq)/Ti 2 . 

Continuity of /+, /_ and their first derivatives at r = r c (8) yields the system of algebraic 
equations for A±, §IB and QB, which has a nontrivial solution if 

J„(+) o 







j ( \ ftn-W go spzjC 1 ) go spetC 1 ) i go cvijC 1 ) 



q+Jn'{+) o -K^i^') ^(fctfS ') 



0. 



(12) 



Here £ = V A2 ( ) " E l = J n(?±r c ), = H^{kr c ), J'(x) = dJ(x)/dx, H'(x) = 

dH(x)/dx. Energies Eq are calculated numerically from Eq. fll2|) . In the first order of the 



perturbation theory, we obtain the energy of the localized states E = E + SE, where the 



energy shift is 



6E —"^i js°ui + r-)<4r ■ (13) 



Density of localized states 

N{E) = 1- T S(E-E U (9)), (14) 

Z7T JO y 

is calculated approximately for E > 0, [ED[ taking 

E u (0) « S o „(0)| cos(20)| + 5£„(0) cos 2 (2#). (15) 



3 Discussion 

Results are illustrated in Figs. 1-3. For a <i-wave "pancake" vortex, only the lowest localized 
state [y = 1/2) is significantly separated from the gap edge for kp^o(0) ~ 1 (Fig. 1). The 
energy shift, Eq. flT3|), of higher states [y > 5/2) is practically negligable. Density of the 
lowest localized states (v < 5/2), calculated within this model, is presented in Fig. 2. Only 
the lowest peak is considerably separated from the gap edge, because of the small kp^o(0). 
One large peak is distinct from the gap edge for 0.7 < k F ^ (0) < 1.2, which may be the 
case for YBCO, and none (gap structure) for fcp£o(0) < 0.7, as for BSCCO, which is in 
consistency with larger gap for the latter (Fig. 3). This can be the simplest explanation of 



the experimental observations. [12, 13 



In conclusion, an approximate analytical method is developed in the same manner 
as the Caroli, de Gennes and Matricon approach, adjusted for the case of <i-wave pairing 



and small radius pancake vortices, where the quasiparticle spectrum is not quasicontinuous. 
Two-dimensional and strong coupling nature of superconductivity in cuprates makes the 
density of localized quasiparticle states in the vortex core completely unlike the normal- 
metal DOS. 
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4 Figure Captions 

Fig. 1. Quasiparticle energy -Z?/A vs angle 9, in the vortex core of radius r c = £,o(9), for 

^f£o(0) = 1-5. Dashed curves: E Q /A Q . Shaded region: unbound states. 

Fig. 2. Density of localized states in the vortex core of radius r c = £o(9), for k F £ (0) = 1.5. 

Fig. 3. Quasiparticle energy of two lowest bound states for 9 = as a function of A /Ep = 

2/nk F U0). 



7 



References 

[1] C. Caroli, P. G. de Gennes and J. Matricon, Phys. Lett. 9, 307 (1964). 
[2] F. Gygi and M. Schliiter, Phys. Rev. B 43, 7609 (1991). 

[3] J. Bardeen and M. J. Stephen, Phys. Rev. 140, 1197 (1965).; ibid 187, 556 (1969). 

[4] D. Rainer, J. A. Sauls and D. Waxman, Phys. Rev. B 54, 10094 (1996). 

[5] N. Hayashi et al, Phys. Rev. Lett. 80, 2921 (1998). 

[6] N. Schopohl and K. Maki, Phys. Rev. B 52, 490 (1995). 

[7] M. Ichioka et al, Phys. Rev. B 53, 15316 (1996). 

[8] Y. Morita, M. Kohmoto and K. Maki, Phys. Rev. Lett. 78, 4841 (1997). 

[9] Y. Morita, M. Kohmoto and K. Maki, preprint (1998). 
[10] M. Franz and Z. Tesanovic, Phys. Rev. Lett. 60, 4763 (1998). 
[11] H. F. Hess et al, Phys. Rev. Lett. 62, 214 (1989).; ibid 64, 2711 (1990). 
[12] I.Maggio-Aprile et al, Phys. Rev. Lett. 75, 2754 (1995). 
[13] Ch. Renner et al, Phys. Rev. Lett. 80, 3606 (1998). 
[14] G. E. Volovik, Sov. Phys. JETP 58, 469 (1993). 

[15] K. Krishana et al, Phys. Rev. Lett. 75, 3529 (1995).; Science 277, 83 (1997). 
[16] J. A. Clayhold et al, Phys. Rev. Lett. 53, 8681 (1996). 
[17] T. Nagaoka et al, Phys. Rev. Lett. 80, 3594 (1998). 
[18] A. van Otterlo et al, Phys. Rev. Lett. 75, 3736 (1995). 

[19] M. Tinkham, Introduction to Superconductivity (McGraw-Hill, New York, 1996). 

[20] I. Knezevic and Z. Radovic, in Superconductivity and Related Oxides: Physics and 

Nanoengmeering III (Proc. of SPIE, 3481, 1998). 



8 




Figure 1 




Figure 2 




Figure 3 



